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Abstract 

We compute the energy and momentum of a regular black hole 
of type defined by Mars, Martm-Prats, and Senovilla using the Ein- 
stein and Papapetrou definitions for energy-momentum density. Some 
other definitions of energy-momentum density are shown to give mu- 
tually contradictory and less reasonable results. Results support the 
Cooperstock hypothesis. 

1 Introduction 

An important outstanding issue in the General Theory of Relativity is the 
problem of finding an acceptable definition of energy-momentum localization. 

H. Bondi PP has argued that General Relativity does not permit a non- 
localizible form of energy, so in principle, an acceptable definition of local 
energy- momentum density should exist. 

The first attempts to identify such a density by Einstein himself j2] as 
well as Landau and Lifshitz Papapetrou and Weinberg jS], were all 
defined in terms of non-tensorial, coordinate-dependent pseudotensors. Un- 
fortunately, the physical meaning of the pseudotensors has been and remains 
a serious point of contention. (For example, one train of thought suggests 
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that a pseudotensor approach could conflict with the equivalence principle 
(0)). 

Attempts to find alternative definitions of a local energy-momentum den- 
sity (by M0ller |7j, Komar 0, Penrose [H], and others), however, have not 
been substantially more successful; for example, Bergqvist [TU] showed that 
no two of the many coordinate-independent definitions of mass give the same 
result for the Reissner-Nordstrom and Kerr spacetimes. 

By contrast, Aguirregabiria, Chamorro, and Virbhadra p!T| showed that 
for all spacetimes of Kerr-Schild class, the energy-momentum complexes de- 
fined by Einstein, Landau-Lifshitz, Papapetrou and Weinberg all give the 
same answer. In Virbhadra noticed that the abovementioned energy- 
momentum complexes coincide for a class of solutions more general than 
Kerr-Schild. For example, Rosen and Virbhadra ^3] showed that several 
energy-momentum complexes give the same and a meaningful result for the 
Einstein-Rosen metric describing cylindrical gravitational waves; note that 
this metric is not of Kerr-Schild class. 

Furthermore, Chang, Nester and Chen ^21, showed that these 
energy-momentum complexes can be redefined as the value of a Hamiltonian, 
which has a surface and a boundary term. The pseudotensors of Einstein, 
Landau-Lifshitz, Papapetrou, Weinberg and M0ller can then be defined in 
terms of the boundary term of the Hamiltonian. Thus the pseudotensors 
are in fact quasi-local, and under this formulation the energy-momentum 
complexes no longer confiict with the equivalence principle. 

As a result there has been a renaissance of interest in and use of one or 
more energy-momentum complexes in a wide variety of contexts |^. 

In general the various energy-momentum complexes give different answers 
|14j . [T7] . Chang, Nestor and Chen jT^j pointed out that the physical sig- 
nificance of the boundary term is still unclear, and hence there are a priori 
many choices possible for an appropriate boundary term. 

It is therefore natural to ask if any of the various definitions of quasi-local 
energy-momentum density is preferable to the others. 

In this paper we calculate the energy and momentum distributions of 
a class of regular black holes defined by Mars, Martin-Pratt, and Senovilla 
JHl using several of the complexes mentioned above. We will show that the 
Einstein energy-momentum complex is one of the complexes that gives a 
meaningful answer for the energy and momentum distributions, supporting 
the hypothesis that the Einstein energy- momentum complex is a 

reasonable definition of energy-momentum density. 
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We use the convention that Latin indices take values from to 3 and 
Greek indices take values from 1 to 3, and we take units where G = 1 
and c = 1. The comma and semicolon, respectively, stand for partial and 
covariant derivatives. 



2 The Mars— Martm-Pratt— Senovilla regular 
black hole 

Although the singularity theorem of Hawking and Ellis ^H] forbids the exis- 
tence of regular black holes with matter content satisfying the strong energy 
condition everywhere, it is still possible to construct regular black holes satis- 
fying the so called weak energy condition everywhere (the energy-momentum 
tensor Tab must satisfy TabW"'W^ > where W is any timelike vector at any 
point of spacetime.) A famous example of such a regular black hole was given 
by Bardeen [20]. (Also see Sharifi [21 for the energy computation). In this 
paper we consider another class of regular black holes constructed by Mars, 
Martin-Pratt and Senovilla ^H]- 

Consider the following line element: 

ds^ = e^^^'^^xir)du^ -2e^^^'^dudr-r\de^ + sm\e)d(l)^) (2.1) 

where 

2M(r) , , 

X(r) = 1 ^ (2.2) 

r 

The functions M and /5 of r are chosen so as to satisfy the following three 
constraints: 

1. the spacetime must represent a Schwarzchild black hole, (i.e. 

M (r) = m 13' {r) = for r > tq with tq < 2m 

where m is some positive constant which represents the total mass of 
the black hole, and prime denotes the derivative with respect to r), 

2. the model fulfills the weak energy conditions everywhere, and 

3. the spacetime is regular. 
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In the paper ^H] the authors give two exphcit examples of models that 
fulfill all these conditions (hence models do exist). In this paper we choose to 
leave M and 13 general in all of our computations; the computations are all 
done for an arbitrary choice of functions M and /3 of r that satisfy the above 
criteria. However, for the graphs we will use the following solution from 



M(x) = mx^(10 - 15x + 6x^) (2.3) 

where x, < a; < 1 is the adimensional variable x = r/2m, and the corre- 
sponding function /? can be integrated from 

(3'{x) = 5{l-xyx (2.4) 

Note that the line element 1)2.11) is a modification of the 
Eddington-Finkelstein form of a Schwarzschild-like spacetime 

ds^ = -x{r)du^ + 2dudr + r'^^dO'^ + sin2(^)c/02) (2.5) 

by e^^*^''^ in the "m" direction, since the (unorthodox) transformation 
gets rid of this term. Hence, the energy of the interior of a 
black hole defined by the line element ()2.ip should be modified from the 
energy of a black hole defined by (|2.5|) by the scalar field e^^'-^^. 

Thus from physical considerations we expect the energy of the region 
enclosed in a 2-sphere to be given by the product M(r)e^^^''^ of the mass 
function M(r) and the scalar field e^^^'^\ 

In order to calculate the energy and momentum distributions in the next 
two sections we will need to rewrite ()2.H) in cartesian coordinates. 

Transforming the line element ()2.1|) to cartesian coordinates via 

u = t + r (2.6) 

X = r sin(6') cos(0) (2.7) 

y = r sin(6') sin(0) (2.8) 

z = rcos(^) (2.9) 



yields 



ds'^ = —dx^ — dy^ — dz^ H [xdx + ydy + zdzf' 



+ e 



1 - e^'^^'') 

'f — dz" + - 

2/3(r) " 



de - (1 - e2^('-)x(r))Mt + -ixdx + ydy + zdz)f\ (2.10) 

r J 
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Note that (j2.ip (equivalently (|2.1(J|) ) is a static spherically symmetric 
metric, but is not a Kerr-Schild class metric. In fact, (j2.H) falls outside 
the class of metrics for which the energy-momentum complexes of Einstein, 
Landau-Lifshitz, Papapetrou, and Weinberg coincide. We will demonstrate 
this explicitly in the following two sections. 

3 Energy-Momentum using the Einstein Energy- 
Momentum complex 

The Einstein energy-momentum complex is 

= ^H^, (3.1) 



where 

'-9 



eg and e° denote energy and momentum density respectively. The energy 
and momentum components are given by 

P^ = J j j Q^^dx^dx^dx^ = -1- J j H^'^riadS (3.3) 

where we applied Gauss's theorem to get the second equality (Pq and Pa stand 
for the energy and momentum components respectively, and ria denotes the 
outward unit normal to the infinitesimal surface element dS). 

Using ()2.10p in ()3.2p . it is easy to compute the twelve superpotentials 
as follows: 

= =^e'^^^\-l + Xir)) (3.4) 
H'o' = ^e^'^^H-l + Xir)) (3.5) 
H^' = ^e'^'^^H-l + X{r)) (3.6) 

i7?i = l(x2(2-2e2^Mx(r)) 

- {y^ + z''){l - 2r(3'{r) + e^^^'^^x{r){l + 4r/?'(r)) + re^^^'\'{r))) (3.7) 
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Hf = ±(x|/(l-2r/?'(r) + e2^('')x(r)(-l+4r/3'(r))+e'^("Vx'(r 



HI 
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H^^ = l(xy(l-2r/5'(r) + e2^(^')x(r)(-l+4r/5'(r))+e2'^('^Vx'(r; 



iy\2-2e'^'^Mr)) 



- {x^ + z^){-l - 2r/3'(r) + e^^'^'\{r){\ + 4r/3'(r)) + re^^^x'lr; 
= J_(^2(i_2r/?'(r) + e2^('')x(^)(-l+4r/?'(r)) + e2^(''Vx'(r; 

- 2r/5'(r) + e^^^''\{r) (-1 + 4r/3'(r)) + e^f^^'\-y^{r\ 

\ {yz{\ - 2rp'{r) + e^'^^'''>x{r){-l + 4r/3'(r)) + e^^^^'^ rx' [r] 



h: 
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) (3.8) 
) (3.9) 
(3.10) 



(3.11) 
(3.12) 

(3.13) 

(3.14) 



- {x^ + y^){l - 2rl3\r) + e^<''^'\{r){l + Ar(3\r)) + re2^(")x'(r))) (3.15) 

(here prime denotes derivative with respect to r). 
Now applying ()3.3p one gets 

E = e^^^''^M{r) (3.16) 

for the energy component. Note that ()3.16p is the value for E we deduced in 
section 2. Furthermore note that for r > 2m, ()3.16p reduces to E = m, the 
energy expression for a Schwarzchild black hole. 

Using ()3.16|1 and recalling our specific choice of model in section 2, in 
Figure 1 we plot E/m (on the y-axis) versus r/2m. 

Similarly, the momentum components are: 



Pi = P2 = Pg = 

as expected for a static spherically symmetric black hole. 



(3.17) 



6 



0.2 0.4 0.6 0.8 1 
Figure 1: E/m vs. r/2m for ()3.16|1 

4 Energy computations involving other Energy- 
Momentum complexes 

The Landau-Lifshitz energy-momentum complex is 

= Y^A^'^'"^,,^ (4.1) 

where 

X^^"^ = -^(^'V""-//") (4.2) 

and denote energy and momentum density respectively. The energy 
and momentum components are given by 

V'^dx^dx^dx'' = —/"/" A'°"" ^mriadS (4.3) 

lOTT J J 

where we applied Gauss's theorem to get the second equality (Pq and Pa stand 
for the energy and momentum components respectively, and na denotes the 
outward unit normal to the infinitesimal surface element dS). 

Using ()2.10|) in ()4.2|) . it is easy to compute the twelve superpotentials A*°" 
as follows: 

Aooi = ^(i_2e2'3M + e^^Mx(r)) (4.4) 
A002 ^ ^(i_2e2/5M + e^^Mx(r)) (4.5) 
A°°3 = ^(l-2e2^('-) + e^^('-)x(r)) (4.6) 




+ (y' + - 2r/?'(r) + e2/^Wx(r)(l + 4r/3'(r)) + re^^^'Wr)))) (4.7) 

^102 ^ _ J_ (^2.^g2/3(r) ^ ^ _ 2^^/(^) + e2^(")x(r) (- 1 + 4r/3'(r) ) + e'^^'^ Vx'(r) ) ) 

(4.8) 

Ai°=^ = -l(xze'^('')(l-2r/3'(r) + e2^('-)x(r)(-l + 4r/3'(r)) + e2^('-Vx'(r))) 

(4.9) 

^201 ^ _J_(^^g2/3(r)(^_2^^/(^)+e2/3(r)^(^)(_;L + 4r/?'(r)) + e2^Wrx'(r))) 

(4.10) 

A202 ^ ^(e2/3M^2|/2(-l + e2'^Wx(r)) 

+ {x^ + ^2)(_i _ 2rp'{r) + e2/^('-)x(r)(l + 4r/5'(r)) + re^'^^^'Wr)))) (4.11) 

A203^ -l(y^e2^«(l-2r/3'(r)+e2'^('-)x(r)(-l + 4r/3'(r))+e2'^('-Vx'(r))) 

(4.12) 

^301 ^ _J_(^^g2^(r)(i_2^^'(^)+e2/3(r)^^^^)^_;^^4^^.^^^)^^g2/3(r)^^/(^^ 

(4.13) 

^302 ^ _ J_ (^^^g2/3(r) ( ^ _ 2^^/ (^-) + g2/3(r) ^ -) ( _ ;L + 4r/?' (r ) ) + e^/^^'' Vx' (r ) ) ) 

(4.14) 

A303 ^ ^(e2/3W(2^2(_;^ ^ e2^«x(r)) 

+ {x^ + 2/')(-l - 2r/3'(r) + e2^('-)x(r)(l + 4r/5'(r)) + re^^^'^^x'ir)))) (4.15) 



Now applying (j4.3p one gets 

E = -^(1 - 2e2^M + e^^M(l - ^^)) (4.16) 

for the energy component, and 

Pi = P2 = P3 = (4.17) 

for the momentum components . Ahhough formula ()4.1(i|l reduces to the 
correct expression {E = m) for the region r > 2m, the energy component 
fl4.16|) differs from ()3.16|) for r < 2m. Thus fj2.1|) falls outside the class 
of metrics for which the energy-momentum complexes of Einstein, Landau- 
Lifshitz, Papapetrou, and Weinberg coincide. Furthermore, negative energy 
is not physically meaningful for electrically neutral nonrotating massive ob- 
jects, while the energy component (j4.16p is negative for small r. (see Figure 
2, where we plot E/m (on the y-axis) versus r/2m.) Therefore, based on 
physical considerations, ()4.16|) appears to be incorrect for r < 2m, which 
supports the contention that the Einstein energy-momentum complex is a 
more reasonable computational tool. 




0.6 0.8 
Figure 2: E/m vs. r/2m for ()4.16|) 

The M0ller energy-momentum complex is 

J- = l^J-i (4.18) 

where 

J^'' = ^^{9^a,b-9^b,a)9''9''' (4.19) 
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Jq and J° denote energy and momentum density respectively. The energy 
and momentum components are given by 

p. = J j j J^dx^dx^dx^ = 1- j J J^^^'n^dS (4.20) 

where we apphed Gauss's theorem to get the second equahty (Pq and Pa stand 
for the energy and momentum components respectively, and denotes the 
outward unit normal to the infinitesimal surface element dS). 

Using ()2.10|) in ()4.19p . it is easy to compute the twelve superpotentials 
JP" as follows: 

rp2/9(r-) 

X' = -^{'^X{r)f3'{r) + x'{r)) (4.21) 
X' = ^^(4x(r)/3'(r) + x'(r)) (4.22) 
X' = -^{Mr)P'{r) + x'{r)) (4.23) 



J-oi = l(rx2((-2+4e2'3Mx(r))/3'(r)+e2^('-)x'(r)) + (y'+^')(-l+e'''(^)x(r))) 

(4.24) 



^ ^(i_2r/5'(r) + e2'3«x(r)( 

J^' = ^(l-2r/3'(r) + e2^«x(r)( 
J2' = -,{ryM-2+4e'^^^\{r))P'{ry 



-1 + Ar(3'{r)) + re^^^^'^x'ir)) (4.25) 
-1 + 4r/?'(r)) + re^^^''^x'{r)) (4.26) 
-1 + 4r/?'(r)) + re^'^^'^^x'ir)) (4.27) 

e'''(^V'(r)) + (x2+^')(-l+e2^(^)x(r))) 

(4.28) 
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Jf' = ^(1 - 2r/3'(r) + e^<^^'\[r){-\ + 4r/?'(r)) + re'^^^^x'l^)) (4-29) 
Jl^ = ^(1 - 2r/3'(r) + e2^(")x(r)(-l + 4r/3'(r)) + re'^^'^^x'lr)) (4.30) 
J3O2 = ^(1 _ 2r/3'(r) + ^<'^'\{r){-\ + 4r/5'(r)) + r^<'^'^^{r)) (4.31) 

^3°' = ^K((-2+4e^^('')x(r))/3'(r)+e2^Mx'(r)) + (x2V)(-l+e'^^'^^x(^^ 

(4.32) 

Now applying ()4.20|) one gets 

E = :!e^«"(4(l - ^),/3'(r) - ^(^^rY -M(r) 
2 r T" 

for the energy component (Here prime denotes derivative with respect to r), 
and 

Pi = P2 = P3 = (4.34) 

for the momentum components. Ahhough again formula ()4.33p reduces to 
the correct expression [E = m) for the region r > 2m, the energy component 
(j4.33p differs from (j3.16j) for r < 2m . (See Figure 3, where we plot E/m 
(on the |/-axis) versus r/2m.) Therefore, based on comparison with our 
expected answer, ()4.33|) is less reasonable for r < 2m, which supports the 
contention that the Einstein energy-momentum complex is a more reasonable 
computational tool. 

The Papapetrou energy-momentum complex is 

S^'' = T^iV*''™,Zm (4.35) 

where 
and 

77*'= = dza^7(l, -1,-1,-1) 
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0.2 0.4 0.6 0.8 1 
Figure 3: E/m vs. r/2m for ()4.33|1 



and denote energy and momentum density respectively. The energy 
and momentum components are given by 

p. = y J J j:'^dx^dx^dx^ = -1- J J N'^p^^n^dS (4.37) 

where we apphed Gauss's theorem to get the second equahty (Pq aiid stand 
for the energy and momentum components respectively, and 71^ denotes the 
outward unit normal to the infinitesimal surface element dS). Note that 
we can apply Gauss' theorem since the metric ()2.1|1 (equivalently ()2.1()j) ) is 
t ime-indep endent . 

Using ()2.10p in ()4.36p . it is easy to compute the only three nonzero su- 
per potentials A^*0" as follows: 

iV00i = Z^e2'^W(-l + xW) (4.38) 
iV002 = Z^e2/^M(-1 + x{r)) (4.39) 
iV003 = Z^e2^M(-l + x(r)) (4.40) 



smce 

jiOa 



AT'"" = for i > 0. (4.41) 
Now applying ()4.37|) one gets the correct expression 

E = e^P^'^Mir) (4.42) 
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for the energy component, and 

Pi = P2 = P3 = (4.43) 

for the momentum components, which agree with the expected values. Thus 
the Papapetrou energy-momentum complex also appears to be a reasonable 
definition for energy-momentum density. 

5 The Cooperstock hypothesis 

The results of this paper support the Cooperstock hypothesis. The Coop- 
erstock hypothesis [221 states that energy is localized to the region where 
the energy-momentum tensor is non-vanishing. This hypothesis would imply 
that there is no energy-momentum contribution from "vacuum" regions of 
space-time. If true, this hypothesis would have broad implications. For ex- 
ample, the hypothesis suggests that gravitational waves have no energy and 
that current attempts to detect these waves using bar detectors are doomed 
to failure^. 

6 Conclusions 

In this paper we calculated the energy of a regular black hole of type de- 
fined by Mars, Martin-Prats and Senovilla (a "MMaS-class" black hole). We 
showed that its energy and momentum could be computed by using the Ein- 
stein and Papapetrou energy-momentum complexes, and that these calcula- 
tions agree with physical intuition. Furthermore, we showed that the energy 
and momentum of a MMaS-class black hole could not be computed in gen- 
eral using the energy-momentum complex due to Landau and Lifshitz, that 
the energy-momentum complex due to M0ller gives a less reasonable answer, 
and that these two calculations give results inconsistent with each other and 
with the energy-momentum complexes of Einstein and Papapetrou. 
The results of this paper support the Cooperstock hypothesis. 

would like to thank an anonymous referee for pointing out that, properly speaking, 
only detection schemes that use a transfer of energy, such as bar detection schemes, are 
affected by the hypothesis, and that laser inteferometry (for example) may still be a viable 
detection technique (see [3^) 
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Chang, Nestor and Chen jT2] stress the importance of identifying appro- 
priate quasi-local energy-momentum criteria in order to identify the right 
boundary conditions in the definition of energy-momentum density. The re- 
sults of this paper suggest that the Einstein and Papapetrou energy-momentum 
complexes may define reasonable boundary conditions. However, the situ- 
ation is more complicated in general. For example, the Einstein energy- 
momentum complex is not symmetric in its indices and so cannot be used 
to calculate angular momentum. Furthermore, whereas both the Landau- 
Lifshitz and the Papapetrou complexes (as defined in this paper) are sym- 
metric in their indices, hence capable of defining conservation laws of angular 
momentum, the Landau-Lifshitz complex seems to give more reasonable an- 
swers in some cases than the Papapetrou complex (see for example |24j). 
Additionally, the Papapetrou complex seems much less amenable to calcu- 
lation for time dependent metrics than its competitors. In any case the 
physical meaning of the Einstein boundary term, its apparent dependence 
on cartesian coordinates, and whether or not the Papapetrou boundary term 
is deserving of equal footing with Einstein's, is in urgent need of further 
study. 
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